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We formulate a theory of double perovskite coumpounds such as S^FeReOg and S^FeMoOs 
which have attracted recent attention for their possible uses as spin valves and sources of spin 
polarized electrons. We solve the theory in the dynamical mean field approximation to find the 
magnetic transition temperature T c . We find that T c is determined by a subtle interplay between 
carrier density and the Fe- Mo/Re site energy difference, and that the non-Fe same-sublattice hopping 
acts to reduce T c . Our results suggest that presently existing materials do not optimize T c . 
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Identification of a ferromagnet with high spin polariza- 
tion at room temperature and stable surface properties 
is an important goal in the field of magnetic materials. 
Such a system would allow, for example the fabrication 
of 'spin-valve' devices of greatly improved efficiency for 
magnetic field sensing [ 1 1 , the development of new mag- 
netic recording media [2|, and perhaps the construction 
of improved sources of spin-polarized electrons for 'spin- 
tronic' applications One promising family of mate- 
rials are the 'double perovskites' These are com- 
pounds of chemical formula ABB'Oq, with A an alka- 
line earth such as Sr, Ca or Ba, and B,B' two different 
transition metal ions. Double perovskites in which B is 
Fe and B' is Mo or Re are of particular recent interest 
because they seem || to be metallic ferrimagnets with 
very high magnetic transition temperatures and highly 
spin-polarized conduction bands. However, neither the 
physics nor the materials science of these compounds is 
yet well understood. The magnetic transition tempera- 
ture and whether the ground state is metallic or insulat- 
ing vary as A is changed from Ba to Sr to Ca Q . Mis-site 
(B — B') disorder has a pronounced effect g|] 

In this paper we take a step towards a theoretical un- 
derstanding for these materials. We derive a many-body 
Hamiltonian, using band theory calculations || to fix 
important parameters. We calculate the magnetic tran- 
sition temperature T c , and determine how different ma- 
terial parameters affect it. Our results should provide 
guidance in attempts to design double perovskite materi- 
als with improved properties, and an appropriate starting 
point for calculations of other properties. 

Double perovskites have a crystal structure which gen- 
eralizes the familiar ABO3 perovskite structure by hav- 



ing two B-site ions, which in the ideal structure alter- 
nate in a simple two sublattice pattern. The band theory 
has been determined ||. The conduction bands are de- 
rived from transition metal i?-site t2 g d-orbitals, in agree- 
ment with quantum chemical considerations 0. There 
are six conduction bands per spin direction per unit cell; 
roughly, one triplet arises mainly from the d xy ^ xz ^ yz or- 
bitals on the Fe and the other from the same orbitals on 
Mo/Fe. The occupied bands are fully polarized at T = 0. 

Because the near-fermi-surface bands are derived from 
transition metal d-orbitals, we argue that a simple tight- 
binding parametrization of the band theory is adequate. 
We therefore model the ideal compound as a cubic lat- 
tice of transition metal sites, divided in the usual way 
into two interpenetrating sublattices, which we denote 
as A(Mo/Re) and B(Fc). On each site we include 3 t 2g 
orbitals, d xy ,d xz ,d yz . A good fit to the band structure 
requires both first and second neighbor hoppings, proba- 
bly because the Mo/Re conduction electrons come from 
the id shell, which is spatially extended. If only first and 
second neighbour hoppings are considered, then the ti g 
orbital symmetry implies that a given orbital can mix 
only with orbitals of the same symmetry and only with 
orbitals in the appropriate plane [d xy couples only to d xy 
orbitals in the xy plane). 

The band theory is thus a sum of three two dimen- 
sional tight binding models, which we write as H ba nd — 
H on . stte +H^P+H^P+H^P, with a\ m<T creates an elec- 
tron on site i of orbital state m and spin a. The on-site 
term H on ^ S i te consists of an A — B site energy difference, 
which has a spin-independent term A m (independent of 
orbital m in cubic symmetry) and a spin dependent term 
Jcr. 
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We now fix parameters by computing the density of 
states (shown in Fig. [l]) and comparing to band cal- 
culations. The calculated tight-binding band structure 
displays features that depend sensitively on parameters. 
The parameters which best fit ref. || are shown in Fig. [I], 
with uncertainties of 10 — 15% on all quantities. 
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FIG. 1. Density of states from tight-binding model 
and dynamical mean field theory. The solid line shows the 
tight-binding DOS (t ab = 0.5 eV, t aa = 0.18 eV, t bb = 0.05 
eV, A = 2.35 eV, J = 1.85 eV); DMFT fit-1 is the local 
DOS from DMFT for essentially the same parameter values 
as the tight-binding model (Wab = 0.5 eV,W A A = 0.22eV, 
Wbb = 0.05 eV, J = 2.2 eV and A = 2.0 eV). DMFT 
fit-2 is another fit to the tight-binding model ( Wab = 0.4 
eV,WAA = 0.4eV, Wbb = 0.02 eV, J = 2.5 eV and A = 2.0 
eV ). The first DMFT fit gives a T c that is twice as much 
as the second one. 



For the |f aa | < \t ab \ 
minima and maxima are at i 



case of relevance here the band 
(A - J) - A(t aa + t bb )] t- 

As lo is increased 



^[Ht aa -t bb ) + (A 
through A — J, the density of states has a step up. For 
tbb = 0, the hybridization vanishes along a line in k- 
space which intersects the van- Hove point, causing the 
density of states N(oj) to diverge as N(lo w A — J) ~ 
(lo - w*)- 1 / 2 In [I/O -lj*)} as lo -> lo* = A - J from 
above. A nonzero t bb eliminates the square-root diver- 
gence and moves the van-Hove singularity away from 
lo = A — J. Thus the leading edge of the sharp peak in 
the density of states fixes A — J and the upper and lower 
band edges determine t ab and t aa . Changing t bb from 
Q.OheV to changes the best fit t ab from 0.5eV to OAeV, 
and cause the density of states in the lo ~ A region to be 
dominated by the (In a;)/ ' \fuo~ singularity, making it asym- 
metric (unlike the published band structure calculation 
M). For later use we note that the "i aa -only" bandwidth 
is slightly less than the "i a f,-only" bandwidth. 



The band structure is strongly spin-polarized. In the 
majority spin sector, one band (dominantly B i.e. Fe 
states) is full and the other band is empty; In the minor- 
ity spin sector one band (dominantly A, i.e. Re /Mo) is 
partly filled and the other (B) is empty. In other words, 
the B-site adopts a fully polarized filled-shell configura- 
tion, and motion of A-carriers onto B sites depends on 
the spin, providing the connection between carrier mo- 
tion and spin alignment which favors ferromagnetism. 

One sees from ref. ^ that in the majority-spin sector 
the A-derived states form a much narrower band than do 
the A-derived states in the minority sector. This shows 
that virtual processes which take the B-site from d 5 to 
d state are much less important than those which take 
the B-site from d 5 to d 6 . Further, Hund's rules indicate 
that a d 4 or d 6 state of less than maximal Fe-spin must 
be even less favorable than a maximal spin Fe d A config- 
uration. Therefore, the local physics on the B(Fe) site is 
described by 
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where S is a unit vector representing the direction of the 
spins on the Fe. Here we have chosen the zero of energy 
to be the A-level energy and have restricted to cubic sym- 
metry. Our approximation treats the c? 5 -maximal spin d 6 
energy splitting correctly but is only an approximate rep- 
resentation of the unimportant higher energy states. We 
further argue that the small filling and high degree of spin 
polarization of the A( Mo/Re ) site means that we may 
neglect electron-electron interactions on this sublattice. 
Thus we propose the many body Hamiltonian 
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We solve this model via the dynamical mean field ap- 
proximation (DMFT) [Q. This method has been widely 
applied to models, such as the Hubbard or Kondo-latticc 
model, with one atom per unit cell. In this situation, one 
formulates a local problem specified by the local action 
Sioc — S w g(cj)c^ (lo)c(lo) + Hi n t, where g(to) is a mean 
field function to be determined by a self-consistency con- 
dition. To generalize this structure to the situation of 
present interest we introduce two mean field functions, 
a(w) and b(co) for the A and B sublattices respectively. 
We use boldface to denote tensors depending on spin and 
orbital indices. Our choice of interaction Hamiltonian al- 
lows us to integrate over the fermion fields, yielding 



Sioc = Tr In [a] + Tr ln[b + A + JS ■ a] 



(4) 



Here the trace is over the frequency, spin and orbital in- 
dices, and an average over the orientation of the core 
spin must still be performed. We now write the self- 
consistency conditions. Because an electron on site A 
may hop to either another A-site ion or to another B-site 
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ion but the spin and orbital indices are conserved by the 
band part of H, we require three couplings (Waa con- 
necting A site to A bath; Wab connecting A site to B 
bath and conversely, and W BB connecting B site to B 
bath). We define the local Green functions G a = (a -1 ) 

and Gb = Mb + A + JS ■ l?)^ 1 \ (() means average over 

the partition function constructed from S; oc ) and obtain 



a = uj + n - W 2 AA G a - Wl B G b 
b = uj + f i- Wl B G a - W 2 BB G h 



(5) 
(6) 



These equations are the generalization to the present 
case of the widely-used semicircular density of states 
equations HJhJ . They differ from the equations success- 
fully used to describe CMR manganites 10 ll| in that 
the 'core spin' interaction occurs only on the B-sitc. Note 
that because H hop is diagonal in the orbital index, these 
equations are diagonal as well. However, the partition 
function which determines the spin average involves all 
orbitals. We shall consider only ferrimagnetic states with 
magnetization direction specified by a unit vector rh, for 
which it is convenient to express the spin dependence as 
Oa^(w) = ao(u>) + ai(w)m ■ 3 a p (and similarly for b). 

We fix parameters by comparing the calculated T = 
DMFT density of states to the tight-binding one. Our 
fits are shown in Fig. we are able to reproduce the 
general structure, including the upper and lower band 
edges reasonably well, but there are some differences of 
detail. The DMFT always exhibits a square root sin- 
gularity at ui = J — A. Because this is integrable, and 
indeed is above the chemical potential, it is not impor- 
tant for our subsequent considerations. The DMFT is 
more likely to exhibit a gap in the density of states for u> 
slightly smaller than J— A. In the band theory, the A—B 
hybridization vanishes along lines in /c-space ( for H'^° p , 
this is along cos(p x ) + cos(p y ) — ); this feature is ab- 
sent in DMFT. Finally, van Hove singularities are present 
in the band theory but not in DMFT. The "W^A-only" 



T = bandwidth is 4W AA , the "W^s-only" bandwidth 
is 4W 'ab ; comparison to the tight-binding model suggests 
W AB ~ 2W AA . 
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FIG. 2. Variation of magnetic T c with the ratio of 
the hops Waa/Wab , as calculated from Eq.( (?|) using the 
"best-fit" value (J - A) /Wab = 0.4; n = 1. Band theory 
suggests Waa/Wab ~ 0.6. The values of T c calculated in 
the Kondo limit described in the text are shown as stars. 
The upturn at Waa/Wab > 3 arises because the A-band 
begins to mix with the majority B-states. 

We now evaluate the magnetic transition temperature, 
T c , which we define as the temperature at which the para- 
magnetic solution becomes linearly unstable to a ferri- 
magnetic one. The spin structure of the mean field pa- 
rameters is expressed by the non- magnetic (ao) and the 
magnetic (ai) components, and at T c nonzero values of 
oi, b\ become possible. Linearizing the equations in eti, b\ 
we obtain for T r : 
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Our two different DMFT fits to the tight-binding 
model (Fig. 0), yield T c 's of 495X(DMFT fit 1) and 
200i ; s:(DMFT fit 2), for the band filling appropriate for 
Mo. Eq.( ^) is a mean-field result. Our previous ex- 
perience comparing DMFT to Monte Carlo results for 
double-exchange models shows that the mean-field ex- 
pression overestimates T c by about 25% JlOj] , implying 
physical T c 's of 370 and 150.fr for these two fits. It is in- 
teresting that T c 's are comparable to or larger than T c 's 
observed in the manganese perovskites. We attribute this 
to the higher orbital degeneracy (which allows each and 
to be partially filled, favoring ferromagnetism) and very 



strong Fe-Mo/Re overlap. 

We see that the calculated T c 's depend sensitively on 
parameters. To elucidate this, we study this dependence 
in more detail. We begin with the effects of varying Waa 
shown in Fig. |[ considering first the limit Waa ^ Wab- 
In this case one has a wide band of A site carriers, weakly 
coupled to a spin-polarized "resonant level" on the B- 
site, a situation familiar from heavy- fermion physics [ jl2| . 
Known results for this problem imply a carrier-spin cou- 
pling I - W\ B /{J - A) and a T c ~ I 2 /W A a- The 
DMFT equations can be solved in this limit, leading to 
precisely this behavior (shown as stars in Fig. pi). Thus 
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T c is maximized at Waa — 0. 
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FIG. 3. Variation of magnetic T c with band filling n, in 
the (A + J) — > oo limit, for different values of J — A and 
Waa/Wab = 0.6; Wbb = 0. 

We now turn to the effects of varying the Fe-<i 6 /Mo 
energy difference J — A, and the carrier density, n. 
We measure the carrier density as number added be- 
yond the Fe d 5 configuration; thus Mo corresponds to 
n = 1 and Re to n = 2. Numerical results are shown 
in Fig. |^; a complicated non-monotonic behaviour is ev- 
ident, which may be understood by consideration of the 
simplifying limits Waa — * and (A + J) ^> Wab- 
In this " double-perovskite-double-exchange" limit carrier 
motion is strongly constrained by the need to have car- 
rier and core spins parallel. Now specialize further to 
the limit | J — A| ^ Wab, so that the carriers are es- 
sentially confined to one sublattice, with effective band- 
width W\ B /\ J — A|. The model then maps onto the 
standard double-exchange one, with three orbitals per 
cell and Hund's coupling J + A. The standard argu- 
ments then show that at low n one has a ferromag- 
net with a T c determined by bandwidth (so increasing as 
| J — A| decreases), while (one electron per 

orbital) some other non-ferromagnetic (probably incom- 
mensurate cf p0|) arrangement of spins becomes favor- 
able. The density range over which a non-ferromagnetic 
ground state is preferred increases as J — A decreases, 
essentially because the greater effective bandwidth in- 
creases the energy gain from populating 'wrong-spin' or- 
bitals, and this accounts for the decrease in the ferromag- 
netic T c with decreasing | J — A| at larger n. This trend 
is in qualitative agreement with the variation in T c on 
changing Mo (n = 1) for Rc (n=2) in Ba and Sr-based 
double perovskites p jl^ , [l^] . 

To summarize, we have formulated a many-body 



Hamiltonian which contains the essential physics of the 
Fe-based double perovskite compounds, and have used 
it to determine the factors affecting the ferro(or ferro- 
magnetic transition temperature. Same-sublattice (Mo- 
Mo) hopping is not small in the actual materials, and acts 
to lower T c . The higher band fillings of the Re compounds 
make it more difficult to realize high ferri-magnetic T c s, 
essentially because o fother competing forms of magnetic 
order. We therefore suggest that a search for higher T c 
materials focus on 4c? 1 systems and on ways to reduce 
the same sublattice hopping, as well as more to closely 
match the Fe g! 6 and B' site energies. On the other hand, 
systems based on Ad 2 ions are more likely to exhibit inter- 
esting many-body physics and non-trivial ground states. 

The Hamiltonian and approximation scheme we have 
proposed allows a number of immediate extensions. From 
the tight-binding band-theory and the DMFT self energy, 
the dc and optical conductivities may be calculated. The 
effects of a lattice distortion which lifts the degeneracy 
of the three t^g orbitals can be computed by allowing 
the site energy A to depend on orbital index. Finally, 
the "cavity field" interpretation of dynamical mean field 
theory allows a straightforward investigation of the ef- 
fects of mis-site ( Re-Fe interchange) disorder. Work in 
all of these directions is in progress. 
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